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A major obstacle for the experimental realization of a supersolid phase with cold atomic gases in 
an optical lattice is the weakness of the nearest-neighbor interactions achievable via magnetic dipole- 
dipole interactions. In this letter, we show that using a large filling of atoms within each well the 
characteristic energy scales are strongly enhanced. Within this regime, the system is well described 
by the rotor model, and the qualitative behavior of the phase diagram derives from mean-field 
theory. We find a stable supersolid phase for realistic parameters with chromium atoms. 
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Magnetic dipole-dipole interactions offer a remarkable 
opportunity to explore quantum phenomena with long 
range interactions in cold atomic gases. Of special inter- 
est are atoms with large magnetic dipole moments such as 
Cr, where the influence of the dipole-dipole interactions 
on the atomic cloud as well as the dipole induced collapse 
have been observed [lj-Q- In the presence of an optical 
lattice, the system gives naturally rise to extended Hub- 
bard models with nearest-neighbor interactions (ij, and 
many remarkable quantum states have been predicted 
0-0]. A major obstacle towards the experimental real- 
ization of these states is the weakness of the nearest- 
neighbor interaction due to the magnetic character of 
the dipole-dipole potential, and correspondingly the ex- 
tremely stringent requirements on temperature, trapping 
potentials, and life-time of the atomic system. In this let- 
ter, we propose an experimentally realistic setup for the 
realization of a supersolid phase in cold atomic gases with 
magnetic dipole-dipole interactions. 

A supersolid phase combines two seemingly contradic- 
tion properties, which in most materials appear in com- 
petition with each other: the arrangement of the par- 
ticles in a crystalline structure with a supcrfluid trans- 
port of the particles 0]. While recent experimental ob- 
servation of a superfluid response in solid He are still 
controversial [10l4l2l |. various models in lattice systems 
have extensively been studied in the past and the exis- 
tence of a supersolid phase has been demonstrated using 
quantum Monte Carlo simulations. Of special interest 
is the appearance of a supersolid phase in a triangular 
lattice (13j and the stabilization of a supersolid phase in 
a square lattice by dipole-dipole interaction Q. These 
Hubbard models can be naturally realized with cold po- 
lar molecules or atomic gases with magnetic dipole-dipole 
interactions. While the understanding of the microscopic 
derivation of the Hubbard model for cold atomic gases is 
well understood the nearest-neighbor interactions 
obtained for atoms in a characteristic optical lattice is 
well below a nano Kelvin, i.e., for Cr with strong mag- 
netic dipole moment V ~ 0.2nK with lattice spacing 
a = 500nm - such stringent requirements have not yet 
been achieved in experiments. 



In this letter, we demonstrate that the critical temper- 
ature for the supersolid phase can be increased by several 
orders of magnitude. This opens up a way towards the 
experimental realization of supersolids with cold atomic 
gases using the magnetic dipole-dipole interactions with 
current experimental technologies. The main idea is to 
allow a high filling factor of atoms per lattice site; then 
the influence of the dipole-dipole interaction is enhanced 
by the number of atoms within each well. In the extreme 
situation, the system is then described by a coupled array 
of Bosc-Einstcin condensates. 
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FIG. 1: Phase diagram for C/~ V: The Mott insulating/solid 
phases appear at commensurate fillings (solid bars) and the 
critical hopping at the tip of the incompressible lobes scales 
as t c ~ 1/n. In turn, the transition from the superfluid into 
a supersolid phase behaves as t SB ~ n, i.e., for large fillings 
n ^> 1 the influence of quantum fluctuations is reduced and 
the transition from the superfluid into the supersolid phase is 
well described within mean-field theory. 

For low filling factors, the competition between the in- 
stability towards a supersolid phase and the Mott insu- 
lating phase for integer fillings gives rise to a rich phase 
diagram. Within this regime mean-field theory fails to 
correctly reproduce the phase diagram and it is manda- 
tory to resolve to exact methods such as quantum Monte 
Carlo simulations On the other hand, for high fill- 
ing factors the critical hopping t c for the phase transi- 
tion towards the solid phases/Mott insulators scales as 
t c ~ !/«., while the instability from the superfluid to- 
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wards the supcrsolid phase appears at much higher hop- 
ping with t ss ~ n, sec Fig. [T] Consequently, the com- 
petition between the different phases is reduced and the 
phase diagram is well described within mean-field the- 
ory. In the following, we derive the transition towards 
the supersolid phase and find the optimal parameters for 
the experimental realization in cold atomic gases with 
magnetic dipole-dipole interactions. 

We start with the description of the Hamiltonian for 
the cold atomic gases with magnetic dipole-dipole in- 
teractions. The system is confined into a quasi two- 
dimensional setup with an additional optical lattice 
within the plane. Each lattice site is occupied by many 
particles, and therefore, gives rise to a quasi-condensate 
on each lattice site. Note, that we assume a lattice spac- 
ing larger than the usual spacing in optical lattices in 
order to reduce losses from three-body recombinations. 
For such high filling factors, the validity of the standard 
Hubbard model breaks down. Nevertheless, the system 
exhibits tunneling between the different wells, and an 
interaction term accounting for deviations around the 
mean particle number n. Then, the Hamiltonian is well 
described by the rotor model 
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where Srii = rii—n describes the deviation from the mean 
particle density n within each well, while <pi denotes the 
phase within each well satisfying the commutation rela- 
tion [rii, <j)j] = iSij. The interaction term Vij contains an 
on-site interaction U for i = j, and obeys the charac- 
teristic decay Vi-tj = ya 3 /|Ri — Rj| 3 for dipole-dipole 
interactions with a the lattice spacing and Ri the lat- 
tice vectors. Note, that the rotor model derives from 
the Hubbard model in the limit of large filling factors, 
however, the rotor model remains a proper description 
of bosonic atoms in an optical lattice even in the regime 
where several higher bands are occupied. Its phase di- 
agram has previously been studied close to half filling 
n ~ 1/2 [lj| 17 1 . Here, we first derive the phase diagram 
of the rotor model at large filling n S> 1, and present 
the effective parameters for a realistic experiment with 
chromium atoms in a second step. 

The system is in the superfluid phase for dominant 
hopping with t 3> V,U, and is characterized by a fixed 
phase 4> within each well and a homogeneous particle den- 
sity n. Its mean-field energy per lattice site reduces to 
E /N = (U + Vxo)n 2 /2 - 2tnz with z = 4 the number 
of nearest-neighbors and N the number of lattice sites. 
Here, Xk = J2j^o exp(ikR J -)/|R J | 3 denotes the dipole- 
dipole interaction in momentum space. The transition 
towards the supersolid phase for increasing interactions 
is signaled by an instability in the excitation spectrum. 
Expanding the Hamiltonian to second order in fluctuat- 
ing fields 5n% = Ui — n and 5<f>i — 4>i — (f> around the 



mean-field values, and introducing the momentum repre- 
sentation with 5<j>k = cxp(ikRj)S(f>j /yN and analog 
for Jrik, we obtain the Hamiltonian H a describing the 
excitation spectrum above the superfluid ground state 

H {1 = ^2 ( £ k <50k&/>-k + Vk<5nk<5n_kJ = E k ala k , 



with the single particle dispersion relation £k = 
2tn [z — 2 J2 a cos (ke Q )] and the effective interaction 
^k = e k /(2n) 2 + U (1 + 7 X k) /2. Here, 7 = V/U dc- 
notes the ratio between the strength of the dipole-dipole 
interaction V and the on-site interaction U, whereas e a 
accounts for the unit vectors in direction of the nearest- 
neighbor lattice sites. Introducing the creation (annihi- 
lation) operators a k (ak) for the excitations above the 
superfluid ground state through = i/3k<50k+<^-k/2/3k 
with = ek/4Vk, we obtain the excitation spectrum of 



the superfluid phase E k = (4Vk£k) 
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FIG. 2: a) Shape of the Fourier transformed dipole-dipole 
interaction Xk plotted in the M-direction. b) Superfluid dis- 
persion relation for different values of 7. At 7 = j c the Ro- 
ton minima occurs. The inset shows the different directions 
within the Brillouin zone. 

Next, we analyze the stability of the superfluid phase 
by varying the ratio 7 between the dipole strength and 
the on-site interaction. The result strongly depends on 
the lattice geometry. Here, we focus on a square lat- 
tice with lattice spacing a; the generalization to arbitrary 
lattice structures is straightforward. The quantity \ k is 
maximal at zero momentum with xa ~ 9.02771, while 
it turns negative and minimal at the edge of the Bril- 
louin zone with K = (Tr/a,n/a) and xk ~ —2.64589, 
see Fig. ^l. As a consequence, the excitation spectrum 
exhibits a Roton minima for increasing dipole-dipole in- 
teractions, see Fig. [2b, and eventually turns zero Ek = 
at the critical value 
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Hence, the superfluid phase suffers an instability at 7 C 
via the nucleation of excitations with momenta K. In 
general, one expects that these excitations form a second 
condensate and give rise to a density modulation for the 
system. The novel ground state is characterized by a 
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supcrfluid response due to the condensates at k = and 
K, and a solid order, i.e., the instability signals a phase 
transition from a supcrfluid into a supersolid. 

Now, we analyze the stability and the ground state 
properties of the novel phase by a mean-field ansatz with 
two condensates. For this purpose, we introduce a den- 
sity modulation fij = |c + dcxp [i (KRj + 9)]\ 2 with the 
constraint c 2 + d 2 = n. The density modulation appears 
as the inference of the two condensates at k = and K 
with the relative phase 9, and exhibits a checkerboard 
structure. Inserting this ansatz into the rotor Hamilto- 
nian ([1]), we obtain the energy per lattice site 

= El+4ztd 2 -2d 2 (n - d 2 ) cos 2 9 (|xk| V-U). 

(3) 

Here, Eq is again the mean-field energy of the supcrfluid 
phase. The second term accounts for an increase of ki- 
netic energy due to a reduction of coherence between the 
different wells, while the last term describes the lower- 
ing of the interaction energy via the density modulation. 
The ground state is obtained by minimizing the energy 
with respect to the density modulation and the relative 
phase between the two phases. The optimization of the 
phase requires 9 = 0,ir, which corresponds to the two 
degenerate ground states reflecting the broken (discrete) 
translational symmetry. In the following, this phase will 
be absorbed into the sign of d. On the other hand, the 
condensate fraction d for the mode K exhibits the typi- 
cal Ginzburg Landau behavior for a second order phase 
transition and predicts the lowering of the ground state 
energy for 7 > 7 C , i.e., the appearance of the supersolid 
phase within mean-field theory coincides with the insta- 
bility of the excitation spectrum. The gain in energy via 
the formation of the supersolid phase takes the form 
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A special property of the energy gain is, that it scales 
with the square of the number of particles per lattice site 
n. The supersolid phase is characterized by a checker- 
board density modulation with the order parameter A 
defined by the correlation function 
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for \i — j\ — > 00 and a superfluid density n s describing the 
supcrfluid flow. Within mean-field theory, these quanti- 
ties reduce to n s = 2tz / {\xk\V — U), and A 2 = l-n 2 /n 2 . 
Note, that within a supersolid, the superfluid density is 
reduced compared to the averaged density, i.e., n s < n 
even at zero temperature due to the additional solid 
structure 

The most remarkable result in Eq. Q is the scaling 
of the gain in energy via the formation of the supersolid 
phase with the number of particles in each lattice site, 



i.e., (E ss — Eq)/N ~ n 2 V/2. This energy serves as an 
estimate for the critical temperature T solid for the for- 
mation of the solid structure, i.e., at optimal values the 
critical temperature T solid ~ n 2 V/2 is strongly enhanced 
compared to the single particle nearest-neighbor energy 
V. This observation is in agreement with the recently 
predicted formation of solid structures in very large su- 
perlattices The detailed estimation of the energy 

scales for realistic experimental setup is shown below. 

Finally, we check the stability of the supersolid phase 
against quantum fluctuations and determine the excita- 
tion spectrum above the supersolid ground state. Again 
we introduce the fluctuating field operators 8<f>i = <pi — <$ 
and Srii = m — rij, where n,/n = I + A cos(K Rj) de- 
notes the mean particle density per lattice site within 
the mean-field theory exhibiting a checkerboard struc- 
ture. Then, the Hamiltonian expanded to second order 
in these operators reduces to 
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Note, that the third term involves the modulated density 
with wave vector K. As a consequence, this introduces a 
coupling for modes with momentum k and k + K, which 
becomes obvious in the moment representation 

H ™ = £ S ^kH-k + Ck Sn k Sn_ k + rj Sn k Sn_ {k+K) ] , 



with the parameters £k = 2tn s [z — 2 ^ Q cos (ke Q )], 
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and r\n s = — zt A/(l — A 2 ). This coupling is a result 
of the broken translational symmetry in the supersolid 
phase, which reduces the Brillouin zone. Hence, we ob- 
tain two modes for each momentum value k with the 
dispersion relation, see Fig. [3l 



= 2 (£kCk + Ck+KCk+K) 

± 2 J (£ k Ck " a+KCk+K) 2 + 4£ k £ k+K v 2 - 



(7) 



The lower branch of the dispersion relation accounts for 
the acoustic modes with a linear sound mode for small 
values k, while the second branch accounts for density 
fluctuations of the checkerboard order. Due to the dis- 
crete translational symmetry this mode is lifted to a finite 
value for small values k and corresponds to an optical 
mode. From the dispersion relation Eq. ([7]) , we find that 
the supersolid phase is stable. 

Finally, we can estimate the relevant parameters for 
the experimental realization of a supersolid phase for 
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FIG. 3: Supersolid dispersion relation, given in different di- 
rections of the reduced Brillouin zone, as shown in the inset. 



chromium atoms in an optical lattice [19(. The limit- 
ing energy is given by the nearest-neighbor interaction 
V for magnetic dipolc-dipole interactions. These in- 
teractions are characterized by the length scale add — 
Ho{j, 2 m/ (12irh 2 ) « 0.8nm, where \x is the permanent mag- 
netic dipole moment. Then, the strength V of the dipolc- 
dipole interaction reduces to 
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with E r = ir 2 h 2 /2ma 2 the recoil energy and a the lat- 
tice spacing. Starting with a conventional density for a 
BEC of cold atomic gases, we obtain a filling with n ~ 40 
for a lattice spacing of a ~ Ifim within the plane and a 
prolate shape of each well along the perpendicular direc- 
tion with aspect ration A « 1/4. The nearest- neighbor 
energy reduces to V ~ 0.5Hz, while the characteristic 
temperature scale for the formation of the solid struc- 
ture reaches T SOIid ~ n 2 V/2 m 0AE r . On the other hand, 
the on-site interaction within each well derives as the 
change of energy for the local condensate within the well 
by adding/removing particles, i.e., U = d 2 E loca i[n]. Tun- 
ing the s-wave scattering length a s across the Feshbach 
resonance present for chromium allows us to compensate 
the contribution from the dipolc-dipole interaction, and 
consequently the on-site interaction can be tuned to an 
arbitrary small value. For particle numbers with n ~ 40 
within the quasi-condensate, the influence of the interac- 
tions is weak on the ground state wave function, which is 
dominated by the contribution from the trap and kinetic 
energy. Then, the ground state wave function is well 
described by a Gaussian wave function and the on-site 
interaction reduces to (see Ref. [l9[ for details) 
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where uJ is the mean trap frequency and a^ the cor- 
responding harmonic oscillator length, while /(vA) de- 
notes a dimensionless function with /(1/2) w 0.5. The 
result is derived within first order perturbation theory in 
the small parameter add,n/aho- The stability against col- 
lapse of the quasi-condensate is well guaranteed for such 
small particle numbers. 



The last remaining parameter is the tunneling en- 
ergy t. From the critical value 7 C in Eq. ([2]). we find 
that the tunneling has to be suppressed by the factor 
t < Vn\xK_\/8 rj 6.2Hz. Note, that the allowed energies 
for the hopping term increases again with the number of 
particles. In addition, the superfluid stiffness involves an- 
other factor n. As a consequence, the appearance of the 
superfluid response given by the Kosterlitz-Thouless tem- 
perature T KT ~ H 2 n s /m exhibits the same scaling as the 



transition temperature for the solid, i.e., T K 



l Vj2. 



It is this scaling of the critical temperatures for the solid 
critical temperature as well as the Kosterlitz-Thouless 
temperature, which allows us to improve the experimen- 
tal parameters by increasing the particle numbers within 
each well. The suitable experiments then can be per- 
formed by adiabatically ramping up the optical lattice in 
a BEC of chromium atoms at positive s-wave scattering 
length a s until the proper filling and hopping energies are 
reached. Within this method, t emp eratures well below 
the recoil energy can be reached [20j . Then, lowering the 
s-wave scattering length allows one to pass through the 
supersolid instability and the additional formation of the 
solid structure. 

Discussions with T. Pfau are acknowledged. The work 
was supported by the Deutsche Forschungsgcmeinschaft 
(DFG) within SFB/TRR 21. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 

[7] 



[9] 
[10] 

[11] 
[12] 

[13] 



[14] 
[15] 
[16] 

[17] 
[18] 

[19] 
[20] 



A. Griesmaier, et al., Phys. Rev. Lett. 94, 160401 (2005). 
T. Koch et al., Nat Phys 4, 218 (2008). 

Q. Beaufils et al., Phys. Rev. A 77, 061601 (2008). 

K. Goral, L. Santos, and M. Lewenstein, Phys. Rev. Lett. 

88, 170406 (2002). 

P. Sengupta, et al., Phys. Rev. Lett. 94, 207202 (2005). 

B. Capogrosso-Sansone, et al., Phys. Rev. Lett. 104, 
125301 (2010). 

L. Pollet, et al., Phys. Rev. Lett. 104, 125302 (2010). 
I. Danshita and C. A. R. Sa de Melo, Phys. Rev. Lett. 
103, 225301 (2009). 

A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

E. Kim and M. H. W. Chan, Nature 427, 225 (2004). 
J. D. Reppy, Phys. Rev. Lett. 104, 255301 (2010). 

B. K. Clark and D. M. Ceperley, Phys. Rev. Lett. 96, 
105302 (2006). 

S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 
(2005); D. Heidarian and K. Damle, Phys. Rev. Lett. 95, 

127206 (2005); R. G. Melko et al., Phys. Rev. Lett. 95, 

127207 (2005). 

D. Jaksch, et al., Phys. Rev. Lett. 81, 3108 (1998). 

F. Hebert et al., Phys. Rev. B 65, 014513 (2001). 

A. van Otterlo and K.-H. Wagenblast, Phys. Rev. Lett. 
72, 3598 (1994). 

E. Roddick and D. Stroud, Phys. Rev. B 51, 8672 (1995). 
T. Lahaye, T. Pfau, and L. Santos, Phys. Rev. Lett. 104, 
170404 (2010). 

T. Lahaye, et al., Rep. Prog. Phys. 72, 126401 (2009). 
I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 



